Asteroid Size-Frequency Distribution 

(The ISO Deep Asteroid Survey) 


Final Report (Item 25) 

21 September 2001 


Contractor : TerraSystems, Inc. 

2800 Woodlawn Drive, Suite 264 
Honolulu, Hawaii 96822 

Contact : Dr. Edward F. Tedesco (P.l.) 

(603) 659-5620 (voice) 

(603) 659-2982 (fax) 
etedesco@TerraSvs.com 


ISO Contrac t: 
Project Title : 
Reporting Period : 
Report Date : 


NAS5-97215 

Asteroid Size-Frequency Distribution 
Final Report 
21 September 2001 


This is a pre-print of a paper submitted to the Astronomical Journal 



Table of Contents 


ABSTRACT 3 

1. INTRODUCTION 3 

2. THE ISO DEEP ASTEROID SURVEY - IDAS 4 

3. ISOCAM DATA REDUCTION 6 

3.1 Observation characteristics 6 

3.2 Summary of data reduction 6 

4. ASTEROID IDENTIFICATION 8 

4.1 IDAS asteroids 8 

4.2 Known asteroids associated with IDAS sources 9 

5. SUMMARY 10 

ACKNOWLEDGEMENTS 1 1 

REFERENCES 12 

Tables I 4 

Figures 1 3 


2 


The ISO Deep Asteroid Survey 1 * * 


Edward F. Tedesco 

TerraSystems, Inc., 59 Wednesday Hill Road, Lee, New Hampshire 03824 

E-mail: etedesco@terrasys.com 
and 

Francois-Xavier Desert 

Laboratoire d’Astrophysique, Observatoire de Grenoble BP 53, 

414 rue de la piscine, F-38041 Grenoble Cedex 9 France 

E-mail: Francois-Xavier.Desert@obs.ujf-grenoble.fr 
Keywords: minor planets, asteroids; infrared radiation — solar system 

ABSTRACT 

A total of six deep exposures (using AOT CAM01 with a 6” PFOV) through the 
ISOCAM LW10 filter (IRAS Band 1, i.e., 12 pm) were obtained on a -15 arcminute 
square field centered on the ecliptic plane. Point sources were extracted using the 
technique described by Desert, et al. (1999, A&A 342, 363). Two known asteroids 
appear in these frames and 20 sources moving with velocities appropriate for main belt 
asteroids are present. Most of the asteroids detected have flux densities less than 1 
mJy, i.e., between 150 and 350 times fainter than any of the asteroids observed by 
IRAS (ledesco, et al., 2002a, Astron J., submitted). These data provide the first direct 
measurement of the 12 pm sky-plane density for asteroids on the ecliptic equator. 

The median zodiacal foreground, as measured by ISOCAM during this survey, is 
found to be 22.1 ±1 .5 mJy per pixel, i.e., 26.2 ± 1 .7 MJy/sr. 

The results presented here imply that the actual number of kilometer-sized asteroids 
is significantly greater than previously believed and in reasonable agreement with the 
Statistical Asteroid Model (Tedesco . et al., 2002b, Astron J., to be submitted.). 

1. INTRODUCTION 

Most main belt asteroids are found between 2.2 and 3.4 AU from the Sun and at 
ecliptic latitudes less than 20 degrees. Except for the largest asteroids, the actual 
number above a given size is poorly known. For example, estimates of the number of 
main-belt asteroids with diameters larger than 1 km range from 3x1 0 4 to IxlO 7 
(Farniella and Davis, 1994). 

The asteroid size distribution is important because it provides constraints on models 
of the original size distribution of the planetesimals formed in the inner solar system and 
their subsequent evolution. It is also an important datum in modeling the numerical size 
of the population of near-Earth asteroids and accounting for their evolution from the 
main belt into Earth-orbit-crossing orbits. 


1 Based on observations with the Infrared Space Observatory (ISO), an ESA project with instruments 

funded by ESA Member States (especially the PI countries: France, Germany, the Netherlands, and the 

United Kingdom) with the participation of ISAS and NASA. 
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A typical 1 km main belt asteroid has a V magnitude between 21 and 23. It would be 
a straightforward program to survey, at visible wavelengths, all asteroids in given 
regions of the sky brighter than this limit. However, because, for any given distance, 
visual surveys are biased in favor of discovering higher-albedo asteroids, magnitude 
data alone cannot be used to accurately derive asteroid diameters . This is because the 
absolute brightness of an asteroid depends upon its cross section and albedo, and 
asteroid albedos span a range of at least a factor of 25 (0.02 to >0.5). Moreover, there 
may be systematic trends of albedo with size (Tedesco , 1994). 

Observing thermal emission permits us to obtain an accurate distribution of asteroid 
diameters because, unlike the linear dependence with albedo at visual wavelengths, the 
infrared flux is only weakly dependant on the geometric albedo. For example, on 26 
March 2001 the 100 km main belt asteroid 50 Virginia was at a solar elongation of 110°, 
a typical elongation for space-based infrared observations. Virginia’s visual magnitude 
at this time, given its SIMPS (Tedesco , et al., 2002a) diameter of 99.82 km, would be 
14.9 if its visual geometric albedo were 0.02 and 18.4 if its visual geometric albedo 
were 0.50, a difference of 3.5 mag. The 12.0 pm magnitudes under these same 
conditions, would be 2.16 and 2.46, respectively, or a difference of only 0.3 mag. 
Furthermore, the lower albedo actually results in a slightly higher 12.0 pm brightness 
because in this case the asteroid’s temperature would be higher. Thus, an IR survey is 
slightly biased in favor of discovering lower-albedo asteroids. 

To date there have been three space-based infrared surveys in which asteroids 
have been incidentally observed: The Infrared Astronomical Satellite (IRAS ), the 
Midcourse Space Experiment (MSX, Ml! et al., 1994), and the Infrared Space 
Observatory (ISO) spacecraft, reported on here. For a description of the ISO mission 
see Kessler et al. (1996) and for details on the ISOCAM instrument see Cesarsky et al. 
(1996). 

Results on IRAS asteroids are given in Tedesco et al. (1992 . 2002a) and those on 
asteroids observed by MSX by Tedesco et al. (2002c). 

The IRAS asteroid survey is severely incomplete at low flux levels, i.e., below about 
1 Jy, because IRAS could only detect an asteroid in its survey mode if a known orbit 
was available. Thus, although IRAS observed at infrared wavelengths, it was limited by 
the (albedo-biased) visual surveys in which asteroids are discovered and from which 
data their orbits are calculated. (IRAS knowingly discovered no main-belt asteroids due 
primarily to the poor spatial resolution of its detectors.) 

IRAS and MSX serendipitously observed numerous asteroids in the course of their 
nominal missions. However, due to the way in which their observations were conducted, 
only asteroids with known orbits were identified with the infrared sources these 
spacecraft detected. IRAS observed -95% of the sky and MSX about 10%. Although 
the faintest asteroids detected in these surveys have flux densities of about 150 mJy, 
they are in no way complete to this flux level. The ISO asteroid “survey” (discussed 
below) observed about 0.125 sq deg of sky to a completeness limit of ~0.6 mJy. 

2. THE ISO DEEP ASTEROID SURVEY - IDAS 

The goal of this survey was to cover the maximum area of sky to the faintest flux 
limit possible under the constraints imposed by the zodiacal background and the 
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available observing time. The field was selected to be in the ecliptic plane, near the 
upper limit of the ISO solar elongation constraint {i.e., near 106°), and located west of 
the Sun (to facilitate ground-based follow-up). In addition, the field was chosen to lie far 
from the Galactic plane and to contain no known IRAS sources or bright stars. The 
sensor used was ISO’s AOT CAM01 with a 6-arcsecond PFOV and using the ISOCAM 
LW10 filter (IRAS Band 1, i.e., 12 pm). 

Asteroids move and their flux may vary appreciably on time scales as short as 
minutes. Consequently, the exposure time was chosen to freeze asteroid motion on 
each sub-map, where a sub-map is a three-by-three arcminute area (the size of the 
ISOCAM array) in which each point was observed three times. Each sub-map consisted 
of a 30 second exposure sequence 2 at a fixed position followed by a step of one 
arcminute in ecliptic longitude or latitude where the exposure sequence was repeated 
and ending with another one arcminute step in ecliptic longitude or latitude where the 
exposure sequence was again repeated. The total time spent doing each point in a sub- 
map was 90 seconds. Because the apparent rate of motion for main belt asteroids 
(MBAs), under the observing geometry described above, is between 0 and 60 
arcseconds per hour, the maximum angular distance moved during the time required to 
obtain a sub-map is less than 1 ,5-arcsecond. However, each sub-map was sampled 
three times to create the complete map, and the times between successive sub-maps 
varied from 30 sec to 870 sec. Thus, the maximum distance a main belt asteroid would 
move between samples of a given point in the map is 14.5-arcseconds. 

Figure 1 is a schematic diagram of the map coverage. Each box is one arcminute on 
a side with North up and East to the right. The raster began with the 3’x3’ array located 
in the NE corner of the map, as indicated by the heavy lines around the nine cells in the 
upper right of the figure. One exposure sequence was made at this position and then 
the array was moved one arcminute (cell) west. Seventeen exposure sequences were 
made along a line of constant ecliptic latitude. This brought the array to the end of the 
first row. At this point, it stepped South one arcminute, made an exposure sequence at 
this position, and then made 16 one arcminute steps to the East to complete the 
second row. This process was then repeated until the center of the array had scanned 
17 rows. 

See the movie for an animated version of this figure speeded up by a factor of about 
60. As can be seen from the movie, or the numbers in the figure, each cell around the 
outer one arcminute of the map received no more than three exposure sequences and 
those in the one-arcminute border interior to this region no more than 6. All other cells 
in the map received 9 exposure sequences for a total of 180 seconds each. We refer to 
the region with 9 exposure sequences as the region of complete coverage. 

The 15x15 raster map was obtained in the same way but using 15, instead of 17 
steps. The complete coverage area is 15 sq arcminutes for the maps obtained with the 
1 7x1 7 raster and 1 3 sq arcminutes for those obtained with the 1 5x1 5 raster. 

Two maps as described above were made in June 1996 and another four in June 
1997. A total of 13.64 hours was expended in obtaining the observations presented 


2 A “clean” was performed at the start of each map. This flashed the array to remove the memory of the 
previous observation and required about 240 seconds. Next 25 stabilization frames were taken. The 
actual observations consisted of four 5-sec exposures at each point in the raster, plus an average of ten 
seconds to step to the next raster position. 
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herein. Maps 1, 2, 5, and 6 required 2.44 hours per map, while maps 3 and 4 had 
available 1.94 hours each. A consequence of the decision to keep the exposure per 
map point constant over all six maps was that maps 3 and 4 cover less area. 

The intention was to have maps 3 and 4 (the time for which was granted under a 
Supplemental Observing Request) made at least 36 hours after the end of the previous 
map pair. However, they were scheduled less than 12 hours after completion of the 
previous map and by the time the observing schedule was issued it was too late to 
reschedule them. 

The sample consists of six data sets that are now in the public domain labeled as 
Target Dedicated Time (TDT) numbers: 21103003, 21103004, 57200101, 57200102, 
57200407, and 57200408 (corresponding, respectively, to map numbers 1, 2, 3, 4, 5, 
and 6 in Table 1). Thus, they were taken in pairs during two 24-hour ISO orbits (the first 
3 digits in the TDT) separated by approximately one year, on 15 June 1996 and 10 
June 1997 (i.e., on Julian days 2450249 and 2450609, respectively). 

Figure 2 shows the six images obtained after processing using the technique of 
Desert , et al. (1999), which is further described in Sec. 3. Figure 3 shows all point 
sources with SNR > 3.0 extracted from the ISOCAM maps shown in Figure 2. Squares 
outline the areas sampled nine times. The point size is proportional to the flux density, 
which ranges from 0.3 to 12.2 mJy 

3. ISOCAM DATA REDUCTION 

3.1 Observation characteristics 

A typical dataset consists of 1800 readouts, each with 5.1 seconds of integration, 
through the ISOCAM LW10 filter centered at 12 pm with a bandpass very similar to the 
IRAS 12 pm band. The lens wheel was on the LGe6 position providing a ratio of 6 
arcseconds per detector pixel (which is also close to the FWHM of the Airy pattern of 
ISO). The camera detector consists of 32 by 32 pixels, with one column (number 24, 
disconnected before launch) missing, providing a 3.2 by 3.2 arcmin instantaneous field 
of view. The total survey area was covered by making a raster with ISO at positions on 
a 17 by 17 grid with 60 arcsecond (10 pixel) steps and 60 arcsecond line separation. 
Each position was observed for four readouts, i.e., 20 seconds of integration time. With 
the survey redundancy (a factor 9), the total integration time per sky pixel is about 3 
minutes. The median zodiacal foreground, as measured by ISOCAM during this survey, 
is found to be 22.1 ± 1.5 mJy per pixel, i.e., 26.2 ± 1.7 MJy/sr (the error bar being the 
dispersion among the six surveys of the same area). 

3.2 Summary of data reduction 

The raw data consists of a cube (CISP files) of detector readouts (one every 36 
CAM time units i.e., 5.1 seconds) and an ISO pointing history (IIPH) file. The detailed 
data reduction procedure is described by Desert et al. (1999). Here we give a summary 
along with the specific parameters that were used for the present datasets. First cosmic 
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rays are removed by a time line analysis of each pixel. Long duration glitches are also 
removed and a transient correction applied, using the method described by Coulais and 
Abergel (2000). 

The data timeline is then analyzed with a “triple-beam” linear algorithm that 
basically finds, for each camera pixel, the difference between the signal at one raster 
position and the average of the 2 adjacent position signals. The dispersion of this 
difference for different raster positions indicate the true pixel noise of the 
measurements, because, most of the time, it is uncontaminated by sources. Badly 
behaved pixel values (due to glitches and bad triple-beam X 2 ) are discarded by an 
adapted sigma clipping. We project the difference and dispersion on a final sky map, 
using neighbor pixel approximation, with a 2-arcsecond pixel size. A redundancy 
number is also obtained this way. The projection is done by coadding with an optimal 
weighing and a first order array distortion correction. We used the associations with the 
USNO optical catalog (version A2.0 - Monet et al., 1998) 4 to deduce the offset positions 
(up to 7-arcseconds in both directions) to apply to each dataset map (because of the 
so-called lens filter wheel jitter). 

The final map is then searched for point sources in a selected area where the 
redundancy is two or more. As explained by Desert et al. (1999), we iterate an 
algorithm where a candidate source (found with a top hat wavelet) is fit with a 9- 
arcsecond FWHM two-dimensional Gaussian (for the position and intensity) and the fit 
is removed. This algorithm allows measuring source fluxes near undefined pixels 
without underestimating the flux (as aperture photometry would do). The noise in the 
flux measurement is deduced from the noise map and the Gaussian least-square fitting 
algorithm. The absolute fluxes were deduced using the nominal ISOCAM internal unit to 
mJy conversion factor (/.©., by assuming that the factor has not changed with respect to 
pre-flight expectancy), and by applying a correction factor (1.52) to go from our fitted 
Gaussian beam flux to total point-spread-function integrated flux. 

In Table 2 5 , we give the complete catalog of (527) sources that were detected at the 
> 3 a level in any of the six maps. Column (1) is an identification number; columns (2) 
and (3) the J2000 RA and Dec; column (4) the flux density in band LW10, (5) the one- 
sigma uncertainty in the flux density; (6) the signal-to-noise ratio; (7), (8), and (9) are 
quality flags; (10) the Julian Date of the observation (an average of the, up to, nine 
measurements on each point source that are available); (11) a confusion flag; and (12) 
a code indicating whether the source was in the multiply-sampled region. Columns (13) 
through (19) provide data on sources found within 6-arcseconds of a USNO-A2.0 
catalog ( Monet et al., 1998) visible source (as reduced by CDS-VizieR, http://cdsweb.u- 
strasba.fr/CDS.html) . Column (13) gives the red magnitude from the USNO-A2.0 
catalog; column (14) the name from the USNO-A2.0 catalog; (15) the number of USNO- 
A2.0 sources associated with the ISO source; (16) the number from the USNO-A2.0 
catalog; (17) the distance from the USNO-A2.0 catalog source (95% are within 4- 
arcseconds); and (18) and (19) the distances in RA and Dec, respectively, from the 
USNO-A2.0 catalog source, rounded to the nearest arcsecond. 


4 As reduced by CDS-VizieR, http://cdsweb.u-strasbg.fr/CDS.html 

5 Table 2 is presented in its entirety in the electric edition of the Astronomical Journal. 
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A catalog of 63 inertial sources is given in Table 3. These are sources from Table 2 
that are seen at the same position in at least two maps. An average flux, ranging from 
0.34 to 10.8 mJy, and error is given along with the USNO-A2.0 catalog association. 
Some well-detected sources have no optical counterparts. These are probably external 
galaxies or very slow-moving asteroids. 2 

The complete (i.e., multiply observed) survey area is 225 arcmin' (0.0625 deg' ). 
The densities of stars and galaxies (assuming all non-USNO-A2.0 sources are 
galaxies), respectively, are found to be 0.072 ± 0.017 and 0.045 ± 0.013 arcmin' , for 12 
pm flux densities greater than 0.6 mJy (a value close to the 4a level where a is the 
median flux error in the complete area for one dataset). 

4. ASTEROID IDENTIFICATION 
4.1 IDAS asteroids 

Moving sources were searched for in those areas observed in common within a day 
of each other. Twenty objects were found to move significantly (i.e., by more than 6- 
arcseconds over a period ^ two hours). 

Because asteroids are moving sources, the technique described in Sec. 3.2 for 
obtaining the flux from coadded images underestimates their flux. Thus, we derived a 
rate-of-motion dependent correction factor (FDCor) by offsetting the inertial sources by 
various amounts to simulate their motion and then performing the photometry as 
normally on the coadded map. This resulted in smaller flux values as a function of the 
amount offset, to which we fit a second order polynomial (shown in Figure 7), viz., 

FDCor = 1 .001 - 0.0034*x + 0.01 60*x A 2, (1 ) 

where x = RT*30*Rate; and RT = 17 (for a 17 by 17 raster) or 15 (for a 15 by 15 raster), 
30 is the time per sample, and Rate is the apparent rate of motion in arcsec per revisit 
interval. 

Table 2 contains the uncorrected flux densities and Table 4, which presents the data 
on the 20 sources identified as being asteroids on the basis of forming tracks with two 
or more sightings, gives FDCor and the corrected flux densities, ranging from 0.43 to 
5.7 mJy, for each sighting. 

All of the identified asteroids have SNR > 4. The 1996 field contains four tracks in 
which at least one sighting has a flux density (FD) > 1 mJy, and the 1997 field contains 
six such tracks. Each field contains ten tracks, at least one of which has FD > 0.6 mJy, 
the 4-a completeness limit. 

Combining the results from the two fields gives 5 ± 1 probable asteroids with FD > 1 
mJy and 10 ± 1 with FD > 0.6 mJy. 

Normalizing these results, gives 80 ± 16 asteroids with FD > 1 mJy per sq deg at the 
ecliptic plane, i.e., with diameters greater than about 1.4 km at mid-belt, and 160 ± 16 
with FD > 0.6 mJy per sq deg (diameters greater than about 1.1 km at mid-belt). 

Singletons {i.e., a source detected only once) with flux densities above 0.6 mJy may 
also be present in these fields but they cannot be unambiguously identified using these 
data alone. 
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4.2 Known asteroids associated with IDAS sources 

We associated two asteroids found in the 1996 field with known asteroids. Details 
on these sources are given in Table 5, where column (1) gives, as two rows per 
observation, the IDAS asteroid number assigned in Table 4 in the first row and the 
associated asteroid’s designation in the second, column (2) gives the IDAS source 
number from Table 2, columns (3) and (4) give, respectively, the observed RA and Dec 
in the first row and the predicted RA and Dec, obtained using the Horizons software 
(Giorqini et al., 1996) 6 , in the second row, columns (5) and (6) the observed corrected 
flux density and SNR, respectively, and column (7) the UTC of the observation. 

Of the ten probable asteroids identified in this field, those associated with 1999 AQ23 
and 17971 1999 JZ50 are the brightest. The predicted V-band magnitudes of these two 
asteroids at the time of the ISO observations were 18.2 and 18.1, respectively. This 
means that 80% of the asteroids in the 1996 field have V > 18 and, for those with low 
albedos (0.02), the maximum V is about 25. 

Figure 6 shows the observed ISOCentric positions for sources 3 and 1008 (1999 
AQ23) and sources 4 and 1003 (17971 1999 JZ50) from the 1996 ISOCAM map, 
together with the ISOCentric positions for the two known asteroids. According to the 
Horizons documentation: “The database is updated almost daily with new objects and 
orbit solutions. Comet and asteroid orbits are integrated from initial conditions stored in 
the JPL-maintained DASTCOM database 7 .”. Due to the ephemeral nature of these 
orbital elements, we present, in Table 6, those used in the analysis described here. 

The ISOCAM coordinates, for the mean time of the extracted sources, are plotted in 
Figure 6 as triangles. Because the asteroid position is the mean from detections 
obtained over the ~18 minutes required to map an inertial point on the sky, the 
predicted positions are shown as a series of 19 positions at one-minute intervals 
centered on the mid-time of the local map. The small squares centered on the predicted 
position trails are 6” on a side, the size of an ISOCAM pixel used in this experiment, 
while the figure is ~3’ on a side, the size of the ISOCAM array. 

1999 AQ23 moved a distance equivalent to the size of a pixel during the time 
required to complete the raster scan of its position, while 17971 1999 JZ50 moved 
about twice this distance. For the numbered asteroid, 17971 1999 JZ50, the difference 
between the observed ISOCAM positions and the predicted ephemeris positions are 
2.8” and 0.6”. For the unnumbered asteroid 1999 AQ23, the observed positions lead 8 
the predicted positions by about 11” (or 34 minutes in time). 

As noted in Sec 3.2, the astrometric accuracy of the IDAS positions is better than 4” 
for 95% of the sources associating with USNO-A2.0 sources. The accuracy for moving 
sources is undoubtedly less, but probably not by a factor of two to three. The formal 
accuracy of the asteroids’ predicted positions is less than 1" (based upon output from 


6 The "HORIZONS (http://ssd.jpl.nasa.gov/horizons.html) On-Line Ephemeris System" was created and is 
maintained by the Solar System Dynamics Group, Jet Propulsion Laboratory. 

7 ftp://ssd.jpl.nasa.gov/pub/ssd/Horizons_doc.ps - Version 2.80 June 14, 2000. 

8 The lead-time is the difference between the time the known asteroid is closest to the position of the ISO 
source minus the time of the ISO observation of that source. The track of 1999 AQ23, in ISOCentric 
coordinates, passes less than 1" from the observed positions. 
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the Lowell Observatory Asteroid Ephemeris (ASTEPH) Version 1.5 at 
http://asteroid.lowell.edu/cqi-bin/koehn/asteph run on 18 September 2001). 

In spite of the relatively poor agreement in position for 1999 AQ23 we nevertheless 
present the albedos and diameters for both 1999 AQ23 and 17971 1999 JZ50 under 
the assumption that they are associated with the ISO sources indicated. 

The results are presented in Table 7, where column (1) identifies the asteroid; 
columns (2) through (5) give the absolute visual magnitude, H, heliocentric distance, H 
Dist, geocentric distance, G Dist, and solar phase angle, Phase, from the Horizons 
ephemerides for the mid-time of the ISO observations; column (6) the corrected 
observed mean ISO LW10 band (-IRAS) 12 pm flux density and its uncertainty; 
columns (7) and (8) the computed diameter and geometric albedo using the Standard 
Thermal Model (STM - Lebofskv et al., 1986), D-STM and p H -STM, respectively; and 
columns (9) and (10) the computed diameter and geometric albedo using the Near- 
Earth Asteroid Thermal Model (NEATM - Harris , 1998), D-NEATM and p H -NEATM, 
respectively. 

The results from the photometry are ambiguous. Using the given values for H and 
the infrared fluxes, the derived geometric albedos, from either the STM or NEATM 
thermal model, are either unphysical (-0.9) or implausible (-0.65 to 0.71). 

Physically plausible albedos can be obtained, e.g., by assuming that H for each of 
these asteroids is 1 .0 mag higher than that published, or that H is 0.5 magnitude larger 
and the 12 pm flux density is 50% higher than reported here, etc. (see Table 7). 
Changes of this magnitude for H are not uncommon (see Tedesco. 2002a). And a 50% 
underestimate of the ISO flux density is also quite possible. Furthermore, both of these 
asteroids are located in the inner part of the asteroid belt where, at least for asteroids 
with diameters greater than ~60 km, low albedo asteroids make up less than 10% of the 
population (Gradie and Tedesco, 1982). However, other than ruling out low albedos for 
these asteroids (because for these values of H and an albedo of 0.02 the predicted 
infrared flux densities range from 200 to 440 mJy, two orders of magnitude higher than 
observed) an accurate albedo cannot be determined. 

5. SUMMARY 

There are about 1 60 asteroids per sq deg at the ecliptic plane above the ISOCAM 
LW10 band (i.e., IRAS Band 1, 12 pm) detection threshold of about 0.6 mJy. This 
corresponds to diameters greater than about 1.1 km at mid-belt. For the fields observed 
in this experiment the faintest asteroid source extracted has a flux density of 0.432 ± 
0.085 mJy and SNR = 4.3. 

To put these results in perspective, note that IRAS’ limiting sensitivity was about 150 
mJy at 12 pm whereas most of the asteroids detected by ISO are between 150 and 350 
times fainter. 

The Statistical Asteroid Model (Tedesco et al. 2002b) was run twice on a four sq. 
deg field centered on the ISO field, once for the epoch of osculation of the model’s 
orbital elements (14 Oct 1998) and a second time (yielding the results in parentheses 
below) for the date of the June 1997 ISO observations. 

A total of 1,063 (1,638), of the approximately two million asteroids in the Statistical 
Asteroid Model, were present in this four sq. deg field. Of these, 673 (852) had 
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predicted IRAS 12 pm flux densities greater than 0.6 mJy. Thus, the model gives 180 ± 
20 asteroids per square degree with 12 pm flux densities greater than 0.6 mJy, in 
reasonable agreement with the ISO results reported here. 

The Statistical Asteroid Model result is actually a lower limit because the model does 
not yet include the NEAs or asteroids beyond the Hilda group. More importantly, 
however, is the fact that it terminates abruptly at a diameter of 1 km. If smaller asteroids 
were included some of these would have 12 pm flux densities greater than 0.6 mJy if 
they were close to the Earth. Nevertheless, the ISO data imply that the actual number 
of kilometer-sized asteroids is in reasonable agreement with the Statistical Asteroid 
Model. 

SIRTF’s imagers will have ~5’x5’ FOV arrays with 1.2” pixels and sensitivities of 
-0.015 mJy at 8 pm and -0.37 mJy at 24 pm, about an order of magnitude more 
sensitive than the IDAS limit. The results presented here demonstrate that significant 
numbers of asteroids will be present in virtually all deep exposures taken near the 
ecliptic plane. Using the extreme size-frequency distributions and the IDAS asteroid 
sky-plane density implies that there will be between 3 and 30 asteroids in each limiting 
sensitivity SIRTF image. See Tedesco et al. (2002b) for additional discussion on this 
issue. 

Unfortunately, in spite of the high quality of the ISO data, little can be said regarding 
the diameters (and nothing regarding the albedos) of the unknown asteroids detected. 
This is because their orbits, and hence distances and phase angles, are unknown and 
cannot be reliably computed from the ISO positions alone. Reliable diameter 
determinations for these ISO asteroids will have to await their discovery. And, with V 
mags probably fainter than 22 for most, this is unlikely to happen in the near future. 

In order to fully exploit such space-based infrared data, orbital elements of the 
asteroids must be known. If albedos are to be obtained, then visual wavelength 
observations are required as well. The minimum requirements for obtaining asteroid 
diameters in the absence of supporting ground-based observations is that the space- 
based data must sample the asteroid’s thermal spectrum at a minimum of three 
wavelengths bracketing the peak emission and be taken at appropriate intervals and 
with astrometric accuracies sufficient to allow computation of an approximate orbit. 
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Tables 


Table 1 . Map Field Centers 


Map 

RA J2000 

Dec J2000 

Raster 

Ecliptic 

Longitude 

Ecliptic 

Latitude 





2000 Eq 

uinox 

1 

22:08:38.6 

-08:34:10.4 

17x17 

338.000000 

0.002471 

2 

22:08:42.5 

-08:33:52.1 

17x17 

338.016800 

0.001180 

3 

22:08:38.6 

-08:34:10.4 

15x15 

338.000000 

0.002471 

4 

22:08:41.6 

-08:33:56.4 

15x15 

338.012900 

0.001483 

5 

22:08:38.6 

-08:34:10.4 

17x17 

338.000000 

0.002471 

6 

22:08:42.5 

-08:33:52.1 

17x17 

338.016800 

0.001180 


Table 2. All Point Sources Extracted From the Six IDAS Maps 


This table is available only as a machine-readable table in the on-line electronic edition. 
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Table 3. IDAS Field Inertial Point Sources 


RA Dec 

339 . 502145~” -8 . 506749 


339.527498 -8.568980 


339.532663 j -8.578120 


USNO Name 


0750-21270670 T 79; 


Associations 

8096 

1017; 4008; 8037; 9022 


5; 9069 



|15 339.611103 -8.559631 


16 339.615798 -8.662995 

17 339.625113 -8 . 436006 



:i 

339.648022 

:2 

339.651048 


3 1 339.655479 | 

9.65 



-8 . 541135 
-8 . 605006 
-8 . 567884 
91558 
74583 
51154 
32838 
-8 . 518580 
-8 . 5 15432 
09590 
45261 
-8.4 50058 
42555 
92524 
04695 
64683 
02145 
66331 
-8 . 614536 
79018 
28508 
83767 
63930 
-8 . 7 11600 
86091 
55913 
-8.596113 
-8 . 645489 
32215 
-8.488266 
84950 " 
28036 
-8.534375 
80850 


-8.573803 



61 339.810757 -8.539094 


339.826604 -8.624449 




241; 5041 


1125; 8253 



liVUiHilfciridsBn 


0750-21272547 7; 


1080; 

134; 


13 


55; 


■ifcidtSAfcvjcHiri 


0750-21273204 


0750-21273242 


0750-21273245 


0750-21273317 


1006; 4002; 5010; 8007; 


8044; 9142 
1053; 4080; 5025; 8023 


1011; 4091; 5007; 8015; 


4149; 5126; 8058; 9054 


; 4009; 5019; 8012; 


4131 


1059; 4038; 5076; 9104 


8067 


18.6 

0750-21273599 

19.0 

0750-21273607 

18 . 6 

0750-21273647 


0750-21273857 I 135 


0750-21273912 1 



0750-21274005 


4015; 5008; 8011; 


4004; 5005; 8006; 


1277; 



■.m.ftvmigi 


0750-21274181 


0750-21274288 


0750-21274307 


■.VM.MMfWH.-l 


9014 


4000; 5000; 8000 


4036; 5016; 8032; 



0750-21276278 1019; 


15 




















































































































































































Table 4. IDAS Asteroid Sightings 


Ast 

ID 

RA 

Dec 

SNR 

Q1 

Q2 

Q3 

JD -2450000 

C 

FDCor 

CorFD 

RT 

1 

3 

339.575977 

-8.577812 

25.1 

1 

1 

4 

249.79034 

1 

1.115 

3926 

17 

1 


339.588994 

-8.580972 

17.2 

3 

4 

4 


1 

1.115 

4225 

mm 

2 


339.647069 

-8.463253 

22.0 

0 

4 

4 


u 

1.481 

5580 

ma 

2 

1003 

339.672558 

-8.454937 

24.3 

0 

4 

4 


a 

1.481 

5559 

ma 

3 

5 

339.595198 

-8.558890 

14.8 

K1 

4 

4 


1 


3111 

ma 

3 

ms 

339.596613 

-8.551524 

13.3 

2 

2 

4 


1 


1997 

ma 

4 

■I 

339.754592 

-8.639367 

■aa 

4 

4 

4 


1 


801 

ma 

4 


339.768698 

-8.638212 

4.0 

■1 

HI 

1 


1 


542 

ma 

5 

amci 

339.700352 

-8.499148 

6.3 

4 

4 

4 

249.77819 

1 

1.092 

841 

ma 

5 

1096 

339.712344 

-8.497636 

6.1 

4 

4 

4 


1 

1.092 

597 

ma 

6 

55 

339.651395 

-8.604998 

7.4 

4 

4 

4 

249.78502 

1 

2.836 

2022 

ma 

6 

1184 

339.702285 

-8.602760 

4.1 

4 

4 

IH1 


1 

2.836 

1231 

ma 

7 

85 

339.715214 

-8.523705 

4.5 

mm 

4 

■n 


1 


974 

ma 

7 


339.752778 

-8.528283 

KOI 

4 

4 

4 


1 


908 

ma 

8 

94 

339.630016 

-8.415515 

mm 

0 

4 

i 

249.78104 

0 

1.008 

781 

ma 

8 

1043 

339.631924 

-8.418823 

4.9 

0 

4 


249.88402 

0 

1.008 

EEzHI 

ma 

9 

■tn 

339.634837 

-8.593815 

5.0 

4 

4 

HI 


1 


heishUA 

ma 

9 

nn 

339.656464 

-8.587599 

4.8 

4 

4 

1 


1 

1.336 


■a 

10 

230 

339.647725 

-8.565410 

5.9 

4 

4 

4 

249.78450 

1 

1.137 

636 

ma 

10 

1090 

339.660784 

-8.558875 

5.3 

4 

■D 

4 


1 

Hkm 

690 

■■ 

11 


339.663681 

-8.665462 

16.4 

1 

1 

4 


1 


2828 

■a 

11 


339.680263 

-8.658790 

24.6 

0 

1 

4 

609.60102 

1 

1.266 

4407 

■a 

11 

■ililill 

339.735623 

-8.635671 

29.5 

0 

1 

4 

609.88002 

1 


5221 

ma 

11 

9013 

339.754758 

-8.626909 

14.6 

0 

4 

4 

609.98160 

1 

1.295 

2791 

ma 

12 

4018 

339.647658 1 

-8.645995 

10.2 

3 

HI 

4 

609.54190 

1 


1141 

ma 

12 

5028 

339.659803 

Emma 

7.8 

0 

HI 

3 

609.60146 

1 


984 

ma 

12 

8014 

339.701216 

-8.625838 

15.0 

■9 

2 

4 


1 

1.159 

1659 

ma 

12 

9005 

339.715405 

-8.619379 

16.5 

4 

4 

4 

609.98432 

1 

1.159 

2321 

ma 

13 

4022 

339.672085 

-8.528938 

7.4 

■a 

HI 

4 

609.53713 

1 

1.047 

778 

■a 

13 

5075 

339.679506 

-8.526589 


4 

4 

3 

609.59736 

1 

■I *Vrm 

961 

■a 

13 

8016 

339.699610 

-8.517709 

9.7 

4 

■a 

4 


■■ 


1154 

■a 

13 

9025 

339.706524 

-8.514754 

8.5 

■1 

HI 

4 

609.98147 

a 


987 

ma 

14 


339.761948 

-8.551948 

KOI 

4 

4 

4 

609.53074 

1 

1.092 

866 

■a 

14 

l&IEI 

339.770892 

-8.545951 

10.0 

4 

4 

4 

609.59120 

1 


1137 

■a 

14 

9023 

339.804980 

-8.513935 

■sa 

HI 

4 

4 


0 

■EM 

1209 

■a 

15 

8161 

339.809759 

-8.565357 

wsm 

■a 

■l 

1 

609.87345 

0 

1.234 

690 

ma 

15 

9064 

339.828591 

-8.565200 

4.7 

Bl 

4 

0 

609.97520 

0 

1.234 

713 

ma 

■a 


339.695200 

-8.541871 


HI 

HI 

3 

609.53601 

u 

1.106 

652 

■a 

mrm 


339.706170 

-8.538451 

mm 

4 

4 

■ 

609.59578 

n 


485 

ma 

Bi 


339.741951 

-8.523027 

mm 

4 

4 

■El 

609.87615 

a 

1.115 

826 

ma 

m 

9045 

339.754062 

-8.517344 

4.8 

4 

4 

i 

609.97850 

n 

1.115 

730 

ma 

17 

4119 

339.630591 

-8.616134 

mm 

4 

4 

i 

609.54245 

a 


1609 

■a 

17 

5219 

339.655302 

-8.567826 

■s 

4 

4 

i 

609.59977 

n 

3.715 

1523 

■a 

18 


339.723015 

-8.643114 

6.9 

■a 

4 

3 


n 

1.149 

762 

■a 

mrm 


339.736023 

-8.639538 

7.0 

0 

4 

■l 


a 

1.149 

812 

■a 

■a 

mil 

339.776996 

-8.625246 


■a 

HI 



1 

1.174 

870 

ma 

Sri 

mssa 

339.792203 

-8.619354 

■H 


HI 

LJ 

609.97887 

1 

1.174 

432 

ma 

Be» 

4137 

339.700181 

-8.610512 

■n 

■1 

HI 

HI 


u 

■EuS 

699 

■a 

■El 

5048 

339.721199 

-8.604013 

■oa 

ma 

■Qi 

4 


n 


930 

ma 

■31 

8010 

339.738208 

-8.687407 

14.9 

0 

HI 

HI 

609.88123 



5686 

ma 

B1 

9017 

339.770672 

mmimm 

6.6 

0 

4 

_JL 

609.98219 

o 


2001 

■a 
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Table 5. Known Asteroids in the 15 June 1996 ISO Field 


(IDAS) Asteroid 

ID 

RA 

(deg) 

DEC 

(deg) 

m 

SNR 

UTC 








(1) 

3 

339.57598 

-08.57781 

4032 

25.1 

06 58 

1999 AQ23 


339.57288 

-08.57761 




(1) 

1008 

339.58899 

-08.58097 



09 25 

1999 AQ23 


339.58596 

-08.58042 











(2) 

4 

339.64707 

-08.46325 

5543 

22.0 

06 45 

17971 1999JZ50 


339.64713 

-08.46247 




(2) 

1003 

339.67256 

-08.45494 

5522 

24.3 

09 11 

17971 1999JZ50 


339.67267 

-08.45481 





Table 6. Horizons Initial Heliocentric Osculating Elements for Ecliptic and 
Mean Equinox of J2000.0 and Epoch = 2001 -Oct-1 8.0000000 (TDB), for 1999 AQ23 
and 2001 -Apr-01 .0000000 (TDB), for 1 7971 1 999 JZ50 1 


Asteroid 

Tp 

q 1 

e 

AP 


i 

H 

1999 AQ23 

2451507.3543131 

2.304351999 

0.1217025 

97.89149 

135.42849 

14.72242 

13.5 

17971 1999 JZ50 

2451472.0746345 

1.887653106 

0.168258407 

167.6773526 

128.5014528 

2.9364673 

14.8 


Vp = Time of perihelion passage (Julian date), q = Perihelion distance (AU). e = Eccentricity, AP = Argument of perihelion (deg), 
Omega = Longitude of the ascending node (deg), i = Inclination (deg), H = Absolute visual magnitude (G, the slope parameter, is 
assumed to be 0.15 for both asteroids). These data were obtained from Horizons on 9 Sep 2001. 


Table 7. Aspect Data, Flux, and Derived Diameters and Albedos for Known 
Asteroids Potentially Associated With IDAS Sources 


Asteroid 

V 

mag 

H Dist 
AU 

G Dist 
AU 

Phase 

deg 

Flux 

mJy 

D-STM 

km 

Ph-STM 

D-NEATM 

km 

Ph-NEATM 

1999 AQ23 

18.2 

2.538 

2.053 

22.57 

4.08 ±0.15 

2.80 ± 0.02 

0.90 ±0.01 

3.30 ± 0.02 

0.65 ± 0.01 

17971 1999 JZ50 

18.1 

1.892 

1.333 

31.01 

5.57 ±0.01 

1.50 ±0.01 

0.93 ±0.01 

1.73 ±0.01 

0.71 ±0.01 

1 Below are results assuming the actual H is one full magnitude higher 


1999 AQ23 

19.2 

2.538 

2.053 

22.57 

4.08 ±0.15 

2.41 ± 0.02 

0.48 ± 0.01 

2.95 ± 0.02 

0.32 ±0.01 

17971 1999 JZ50 

19.1 

1.892 

1.333 

31.01 

5.57 ±0.01 

1.30 ±0.01 

0.50 ± 0.01 

1.54 ±0.01 

0.36 ± 0.01 

t Below are results assuming the actual H is 0.5 magnitude higher and the Flux a factor c 

)f 1 .5 greater. 


1999 AQ23 

18.7 

2.538 

2.053 

22.57 

6.1 2 ± 0.15 

3.00 + 0.02 

0.51 ± 0.01 

3.60 ±0.02 

0.34 ± 0.01 

17971 1999 JZ50 

18.6 

1.892 

1.333 

31.01 

8.36 ±0.01 

1.60 ±0.01 

0.52 ± 0.01 

1.90 ±0.01 

0.37 ±0.01 
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Figure 1. Sample Region Coverage of a 17 x 17 Raster Map. Each 1' x 1' cell shows 
the total number of 20 second integration sets made on that area. 
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a 1996 June 15 05:37:57- 08:02:09 UT 



b 



■, 'r< 3,lv if i V’ - 7 ■ 

1996 June 15 08:03:21 - 10:27:33 UT 





-v ... , ■ 

V fS u ~ At ■*’ V SS J * ■' 1 

c 1997 June 10 00:19:15-02:13:35 
UT 


►. r 5f ,■ ^ 



d 1997 June 10 02:14:39- 04:08:59 UT 




- -J. V ■ - ■ P ■ -H 

1# jft . ri_ * ppv . - ■ _■ :*■ “-4 \ 
e 1997 June 10 09:14:03- 11:38:14 UT 



* ■ . .■ T> -i -■ M ■. ' 

f 1997 June 10 11:39:26- 14:03:38 UT 


Figure 2. ISOCAM maps in ecliptic coordinates: 1996 June 15 (a, b) and 1997 June 10 (c, d e, f). North ecliptic latitude is up 
and east ecliptic longitude is to the right. Maps a, b, d, and fare 17x17 rasters while c and d are 15x15 rasters. 
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337.9 338.0 338.1 

Eciptic Longitude (deg) 

a 1996 June 15 05:37:57-08:02:09 UT 


0.1 - 



- 0.1 - 
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337.9 338.0 338.1 

Eciptic Longitude (deg) 

b 1996 June 15 08:03:21 - 10:27:33 UT 



337.9 338.0 338.1 

Eciptic Longitude (deg) 

c 1997 June 10 00:19:15-02:13:35 UT 



337.9 338.0 338.1 

Eciptic Longitude (deg) 


337.9 338.0 338.1 

Eciptic Longitude (deg) 


337.9 338.0 338.1 

Eciptic Longitude (deg) 


d 1997 June 10 02:14:39 - 04:08:59 UT e 1997 June 10 09:14:03 - 1 1:38:14 UT f 1997 June 10 11:39:26 - 14:03:38 UT 

Figure 3. Point sources from the ISOCAM maps shown in Figure 2. Squares outline the areas sampled nine times. The point 
size is proportional to the flux density, which ranges from 0.3 to 10.8 mJy. 
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337.9 338.0 338.1 

Eciptic Longitude (deg) 

Figure 4. IDAS Inertial Point Sources. Squares outline the areas sampled nine times 
per map. The inertial sky within the union of the two large squares was sampled 36 
times (total integration time 720 sec per map pixel) and that within the two small 
squares 54 times (total integration time 1,080 sec per map pixel). The point size is 
proportional to the flux density, which ranges from 0.3 to 10.8 mJy. 
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337.9 338.0 338.1 

Ecliptic Longitude (deg) 


Figure 5. IDAS Non-inertial Point Sources, i.e., all sources from Fig 3 not plotted in Fig 
4. Sources from the 1 996 field are plotted in the top panel and those from the 1 997 field 
in the bottom panel. Squares outline the areas sampled nine times per map. The point 
size is proportional to the flux density, which ranges from 0.4 to 4.3 mJy. Circles 
indicate sources extracted from Maps 1 (top) or 3 (bottom), triangles those from Maps 2 
(top) or 4 (bottom), squares those from Map 5, and diamonds those from Map 6. 
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1999 AQ23 


339.57 


339.58 


339.59 


R.A. (ISOCentric J2000) 


17971 1999 J250 


o -8.45 


339.65 339.66 339.67 339.68 

R.A. (ISOCentric J2000) 

Figure 6. Observed (circles) versus predicted (triangles) positions for two bright known asteroids in the 
1996 ISO map. Each graph is 3’ on a side, the approximate size of the ISOCAM array. The small squares 
centered on the observed position are 6” on a side, the approximate size of the ISOCAM PFOV used. 
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